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Abstract 

Partons are defined as the quanta in a Fock space description of a 
hadron. Gluon saturation is described in the Weizsacker-Wilhams ap- 
proximation for a large nucleus. The elements of DGLAP and BFKL 
evolution are given with the BFKL equation derived in a large-A^c 
dipole formalism. A more general discussion of saturation is given in 
terms of a dipole scattering on a nucleon or nucleus. Possible evidence 
for saturation at HERA is discussed. 

1 Introduction: Partons, classical fields and a 
simple picture of saturation 

In this section partons are defined as the quanta in the wavefunction of a 
hadron, or dressed quark, when quantization is done on the hght-cone and 



"'^Lectures given at International Summer School "Particle Production Spanning MeV 
and TeV Energies," Nijmegen (Aug.8-20, 1999) and XVII Autumn School "QCD: Pertur- 
bative or Nonperturbative?," Lisbon (Sept.29-Oct. 4, 1999) 

^This work is sponsored in part by the Department of Energy, Grant DE-FG02-94-ER- 
40819 



1 



when light-cone gauge is used. The relationship between this "quantum" 
description and a picture where gluons are viewed in terms of a classical, or 
semiclassical, field in a hadron is then given. Although the more classical 
description has limitations it also provides a useful and intuitive way to 
picture many aspects of a high density parton system in terms of strong field 
configurations. Gluon saturation, in parton language, becomes a limit on 
the strength which the light-cone potential can reach, in the language of field 
configurations. 

1.1 The field and gluon distribution in a quark 

Consider the wavefunction of a quark at order g in QCD. The formula 



is illustrated in Fig.l and represents a dressed quark state in terms of a bare 
quark, the first term on the right-hand side of (1), and a bare quark plus a 
gluon. N is a normalization factor, A = ± are the polarizations of the virtual 
gluon, and C—1---8 are the color indices of the gluon. We have suppressed 
the color indices for the quark and instead view as a 3 x 3 matrix in 
quark color space with the quark state being a 3-spinor in color space. The 
integration in (1) is (Pk — dk^(Pk with k± — 4=(/co ± ks). 




X=± c=l 



(1) 




P-k 



k c 



FIG.l 
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Recalling that in light-cone quantization momenta P+ and P are kine- 
matic while P_plays the role of the Hamiltonian, with P_ = P^/2P+ for an 
on-shell zero mass particle, can be written from light-cone perturbation 
theory as 



.3. /^,c/,^ (P' -k,k{X,c)\gJd'^xq(x)-fi,^q(x)A''(x)\p) 
S{p-p)^,{k)^ + (2) 

where d^x = d^xdx- and where the free gluon field can be written as 



Ki^) = E / / "^'^ [e^(A;)a^(A;)e'^-^-^^+^— '^-^+ + h.c] (3) 

with 

[ai{k'), a'^ik)] = SxyScc'S^k - k') (4) 

and where ^^(A; - k') = S{k - U)5{k+ - k'^). 

It is useful to imagine the calculation done in a frame where p+ is large 
and p = in which case 

In the soft gluon approximation (the small-x approximation) k^/p^ « 1, 
and this implies that k^ » {p — k)_. Thus the "energy" of the softest 
particle, the gluon, dominates the energy denominator in (2). Also, in = 
light-cone gauge 

6j(A;) = (6^,6^6^) = (0,^,6^) (6) 

and, again in the soft gluon approximation, one need keep only the term 
from (6) in (2) and (3). Thus, one finds 

2 k ^27r)32A;+ 
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Problem 1(E): Use (2) along with (3)-(6) and U{p — k)^ixU{p) a; 2p+5(^_ 
to derive (7). 

Now define the gluon distribution of the quark by 

xG,{x,Q'')^ jY,d\e{Q^-e)xS{x-kJp+){i;p\ai{k)al{k)\'iPp) (8) 

which immediately gives 

xG,{x, Q^) = Y.j X S{x - k+lp+)Q{Q^ - k^)ipi'{k)rx{k) 

c,A 

= > / x5ix - k+ p+)G(Q^ - k^). (9) 

V 2 2 i (27r)3 2k+ ^ - > ^ > 

Introducing an infrared cutoff, /i^, for the transverse momentum integral and 
using X^c ^ ^ — — ^2N^ and g"^ = Ana one finally obtains 

xGn{x, Q^) = ^en Q^fi^ (10) 

TT 

Equation (8) is a natural and intuitive definition of the gluon distribution as it 
clearly corresponds to the number of gluons per unit of x in the wavefunction 
having transverse momentum less than Q. The free integration of d'^k in (8) 
for ^ < can be viewed as determining that the gluon, and the quark from 
which it came, are point-like and localized down to a transverse size Ax± — 
2/Q. We note that if one takes xG{x, Q^) = 3xGg{x, Q^) = ^^inQ^/jj,^, one 
obtains a result for the proton which is not unreasonable phenomenologically 
for X ~ 10~^ - 10~^ and moderate Q"^ if is taken to be 100 MeV. 

One can associate a "classical" field with the gluon in the quark according 

to 

<(a;) = I d'p'{^,,\A';,{x)\^,) (11) 



or, for i — 1,2, 



^"^^-7(277)3" ^2h\^- ^ ^ 
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It is then natural to define the "classical" field in momentum space to be 

(13) 



Problem 2(M): Take ^ = in (12) and show that 



Problem 3(M-H): Derive (12) using the Feynman graph in Fig. 2 along with 
the light-cone gauge propagator L>^^(A;) = p^f^^x^.-]^ - where rj- 

V — V+ for any vector V^. 



Problem 4(E): Show that Gq, as given in (9), can be obtained from 



Gq{x,Q') 



P 



V 



p-k 



FIG.2 
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The ie prescriptions in problems 2 and 3 may seem a little unusual. We 
shall have more to say about this later on. The fact that the ie's are different 
in the two fc+ denominators in the propagator in problem 3, however, is 
necessary for consistency. If we orient the gluon momentum according to the 
arrow in Fig. 2 then the —ie corresponds to the case that the t] carries the 
index at the tail of the arrow. This uniquely prescribes how the ze's are to be 
given. Finally, problem 4 shows that the gluon distribution can be written 
in terms of the "classical" gauge potential a quantity intimately connected 
with the integral of the field strength as is clear from problem 2. The 
classical field we have defined here for gluons is identical, except for the 
color matrix to that which one naturally defines for photons in quantum 
electrodynamics. 



1.2 The gluon field and distribution in a proton 

A quark can emit, or absorb, a gluon and remain a quark, however, a proton 
cannot. Thus we generalize the idea of a classical gluon field in a proton [|l|], 
and instead of the result given in problem 2 we take 

Ar^^\x) = - gj d%db.^^^^^Qib.-x.)p%b,b^) (14) 

where p is the quark color charge density operator for the proton. It is 
convenient to introduce a matrix notation for the gauge potential where 
= J2c Y ~ 2tr[Y^]- Thus in matrix notation 

Af{x) = -gJ d'b^^^^^Qib^ - x^)p{b, 6_) (15) 

and 



where 



p{k) = J d%db-e-'^-^-+'^+^-p{b). (17) 
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The integration over b in (14) - (17) goes over the volume of the proton with 
b the usual transverse coordinate and 6_ the light-cone analog of the usual z- 
coordinate. We shall only use (15), or (16), for transverse wavelengths much 
smaller than the radius of the proton in which case we may make p more 
precise by writing [|l| 

< mP^'ib') >= ^5{b - mb- - b'_)pib, 6_) (18) 

for the charge-charge correlator in the proton where p{b,b-) is the quark 
number density normalized according to 

J (fbdb^p{b) = Nc. (19) 

Of course (18) simply reflects the fact that for short transverse spatial sep- 
arations only individual quarks are effective while (19) reflects the fact that 
we are working in a valence quark model @], and that we do not consider the 
evolution in the proton which would produce either a quark sea or gluons 
additional to those which will couple to p through (15) and (16). That is we 
shall only allow two gluons to be emitted or absorbed in a proton. Finally, 
when a;_ is sufficiently negative 

-tr [ d% < Af{x + 6, x_)Af (6, x_) = xG{x, Q^) (20) 
vr J 

with identified with Equation (20) shows that we can recover the 

gluon distribution in a nucleon, as given by the independent gluon distribu- 
tions of the valence quarks, in terms of the "classical" gluon field we defined 
in (14). 

Problem 5(E): Derive (20) 

1.3 Non-Abelian Weizsacker- Williams field and gluon 
saturation 

Now we wish to find A^^^'^ for a large nucleus. We suppose that A^^^ is small 
for a nucleon, but for a large nucleus A^lf \ the light -cone gauge potential, 
will be large and non-linearities will have to be taken into account. It is these 
non-linearities which lead to the idea of gluon saturation. Our approach is 



7 



the following. In a nucleus it is well-known that interactions between the 
nucleons are short range and, to a good approximation, these interactions 
are weak enough, because of the diluteness of nuclear matter, that they may 
be neglected altogether in scattering hadrons and real and virtual photons 
on nuclei. If we express the gluon field of a nucleon in a covariant gauge 
then the gluon field will decrease reasonably rapidly outside a nucleon, even 
without having a model of confinement, due to the color neutrality of the 
nucleon. Thus we are going to assume that the "classical" gluon field of a 
nucleus is additive in the nucleons in covariant gauge and then transform 
that result to light-cone gauge where it no longer will be additive, but where 
the partonic interpretation will be manifest . We begin by noting that in a 
covariant, Lorentz, gauge the classical field of a quark becomes 

A^^^\x) = HkH, A[ = A'_ = (21) 

instead of that given in problem 2. Problem 6(E): Show that the same F+j 

as given in Problem 2 can be obtained from F!^^\x) = — Vt A^^'^^\x). 
The Lorentz gauge potential for a nucleon is then given by 

Al''\x,x.) = -gj '^^n[\x-mp{h,x_) (22) 
with the p{h) = Ec f P"^^) 

same as in (15) and (18). The infrared cutoff, 
/i, in (21) and (22) will disappear when physical quantities are evaluated. 
Now we take the Lorentz gauge potential for a large nucleus simply by adding 
the contributions from the individual nucleons as given in (22). In order 
to transform this result to light-cone gauge we use the general formula for 
changing gauges in QCD 

A,{x) = SA^{x)S'^ - i/g{d,S)S-^ (23) 

where is a (space-time dependent) unitary transformation to be deter- 
mined. In order that A^ correspond to a light-cone gauge potential we require 
y4_|_ = in (23) which gives 

-i-S = -WSA'^. (24) 
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Equation (24) is easily solved as 

/X — 
dx'_ (x, x'_ ) } (25) 
-oo 

where the P symbol in (25) indicates that one should order the matrices 
so that A'_,_(x, x'_) comes to the right (left) of A'_,_(a;, x"J) in case x'_ is greater 
(less) than x"_. The choice of writing the path integral from x'_ = — oo to 
x'_ = x_ corresponds to the choice of ie's in the light-cone denominators in 
problems 2 and 3 and in (16). 

Now that S has been determined in terms of the Lorentz gauge potential 
we go back to (23) to get A'f^\ The first term on the right-hand side of (23) 
is not large since it is just a unitary transformation of the Lorentz gauge 
potential. The large term is 

Af^\x,x.)= -'-{d,S)S'' 
9 

leading to|^ 

M'^\x, x_) = - r~ dx'_S{x, x'_)d,A'^{x, x'_)S-\x, x'_) (26) 

where S is given by (25) with A'_^_ given by (22). 

The evaluation of (26) is a bit technical and we just outline the steps 
hereP]. We begin by assuming our "nucleus" is just a thin slice in the z{x-) 
direction. Then, later, we will assemble the various slices to get the total 
result for spherical nucleus. Then for the thin slice 



-9j^dMk'b\f,Wib,xi){^ 

rd%'db' 



S{x, x'_)diA'l^^^ {x, x'_)S ^{x, x_ 
A 



[ d%'db' 6_ 

1 - — 
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(2vr; 



2 p{b',b'_)in^[\x-b'\fi] 



-9'fcde^ I ''-^H\x-bMp\U,b'_)}. 



(27) 



In arriving at (27) we have anticipated our next step and used (18) for the 
quadratic term in p which leads to the pin^ term in (27) with p{li,b'_) the 
quark number density in the nucleus normalized according to / pd^b = NcA. 
For a thin slice of nuclear matter we have kept only quadratic terms and 
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linear terms in p in (27). Now it is straightforward to find, still for a thin 
slice of nuclear matter, 

^ j (fhtr < Af{h + x)Af{h) >= J dx-d%xG{x,Q^)p{b,x.)- 

■ |l-/_'^rf&-^^^p(6,6^)xG(x,g2)| (28) 

with xG{x,Q'^) = Nc^^inQ'^ / jj'^, the gluon distribution in a nucleon and 
where now p{b,b-) is the normal nucleon number density in the nucleus a 
factor ofl/Nc times the p which appears in (27). It is now straightforward to 
go to a thick nucleus by simply exponentiating the second term in brackets 
on the right-hand side of (19). Thus if we call 

N{x^) = -tr I £h < Afih + x)Af{h) > (29) 

then|, i 

N{x')= fd'b^l^^il-e-^'^'/'] (30) 
where the saturation momentum, Qs, is given by 

Ql = ^-^^^R'-b' pxGix, Q\ (31) 

Problem 7(M-H): Derive (27) and (28), from (18), (22), (25) and (26). 
The derivation of (28) from (18) and (27) is not hard but requires good 
bookkeeping skills along with the result of Problem 8. 

Problem 8(H): Show that for small x^ 

j <eUM (6+x-6')/^l -^^1 (&-fe')/^l -^^^'1 (^-^0/^1 -^^^'1 (^-^'+ai)/i|]p(&', h.) 

TT 

= — — x^p(6, 6_)£nl/a;^. 
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1.4 Interpreting saturation in the McLerran-Venugopalan 
model f 



It is tempting to identify N(x^) with xGa{x, Q"^) as was done in (20) for the 
nucleon. This identification is correct if Q'^ = lies far above Q^, but it is 
not correct for below Q'^. What we have actually been calculating is the 
number density of gluons in the nuclear Fock space wavefunction. So long 
as Q'^/Ql is much greater than one this number density agrees with that 
measured in deep inelastic scattering. However, the gluons in the nuclear 
wavefunction have been rearranged in transverse momentum compared to 
what a simple superposition of nucleon gluon densities would suggest and at 
the semiclassical level we are now working the gluon distribution as defined by 
the operator product expansion, and as measured in deep inelastic scattering, 
is still additive. 
Thus 

is the number of gluons per unit transverse momentum in the nuclear wave- 
function while 

dxG r<Px N'i-l 



d%dH 



is the number of gluons per unit of transverse phase space in the wavefunction 
of the nucleus. 



Problem 9(M): Show that the total number of gluons in the nucleon wave- 
function, the integral of (33) over all b and i, is simply the sum of the number 
of gluons in each of the nucleon wavefunctions even though the distribution of 
gluons in i has been modified from a simple additive rule. In particular, this 
means that there is no shadowing in the McLerran-Venugopalan model. 

The result in Problem 9 shows that saturation and nuclear shadowing are 
very different phenomena. Saturation is apparent in|^ 

< Af{b + x)Af{b) >= - e-^«^/^) (34) 
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which follows from (29) and (30) and which says that the light-cone gauge 
potential never becomes larger than 



the "ultimate" value of the potential allowed before gluonic "interactions" 
move gluons from lower transverse momenta to higher transverse momenta. 
We shall come back to discuss saturation in a more general context, and 
especially for quark densities, in the third chapter of these lectures. 



2 DGLAP and BFKL evolution 

Dokshitzer, Gribov, Lipatov, Altarelli and Parisi (DGLAP) 0)01 evolution 
and Balitsky, Fadin, Kuraev and Lipatov (BFKL)P, ^ evolution tell how to 
go from a (given) coarse-grained description of partons in the wavefunction 
to a more fine-grained description. In the case of DGLAP evolution, which 
is equivalent to the usual renormalization group evolution, one passes from a 
coarse-grained description in Ax_l (or in Q) to a more fine-grained description 
in Ax±{oT in Q), while in BFKL evolution one moves from a coarse-grained 
description in At (or in x) to a more fine-grained, shorter time, description. 

2.1 DGLAP evolution in the leading logarithnmic ap- 
proximation 

Recall, from (9) and (10) we found an expression for xGq{x, Q"^). Clearly this 
expression obeys the equation 

Q2^xG,(a:,Q2) = ^. (36) 

If we interpret (35) as telling us how many gluons come from a single quark in 
a given small transverse momentum bin (see (9)) then we can generalize this 
equation easily to the general distribution of gluons in a proton according to 

Q^t|^xG(x, Q') = ^x'qix', Q') + ^ f' ^x'Gix', Q'). (37) 

OU'^ TX Jx X TT Jx X 
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In the first term on the right-hand side of (37) the factor is exactly 
as in (36) while the factor multiplying is the number of quarks (and 
antiquarks) having longitudinal momentum greater than the gluon emitted. 
The second term on the right-hand side of (37) corresponds to emission of 
the "observed" gluon off higher longitudinal momentum gluons and Nc — Ca 
naturally replaces Cp- 

Problem 10(H): Drop the x'q{x', Q"^) term in (37) and solve for xG{x, Q"^) in 
the small-x and large-Q^ limit, in the leading double logarithmic approxima- 
tion. First suppose a is a constant. Then, if you can, solve when a — 2 

with6=ii^. 



2.2 The diple picture of BFKL evolution [10- 12] 
2.2.1 Lowest order 

At lowest order the cross section for heavy onium-heavy onium scattering is 
given in terms of the amplitude illustrated in Fig.3, where the lines P — ki 
and P' — k[ are heavy quarks and the lines ki and k[ are heavy antiquarks. 




FIG.3 

We find it convenient to describe the onium wavefunctions using trans- 
verse coordinate space variables, x^q, x^i, x^x^i, and longitudinal momentum 
fractions, Zi, 1 indicated in Fig.3. The coordinate and 

momentum space wavefunctions are related by 

^l^^'\x,„z,) = / lhLe^kr^-..^io)^k„z,), (38) 
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where and the superscript (0) indicates a lowest order (no 

soft gluons) wavefunction. In terms of = the lowest order cross 

section is 



J d^xoi J dzi¥^\xoi, zi) J d'^x'oi J dz[^^^\x'oi, z[)add{xoi, x'^^), 

(39) 

where add is the dipole-dipole scattering cross section 

add{x, x') = 27ra^x'{ 1 1 + in— ] (40) 



X 



with x> the greater of x and x' and with a;< the lesser. Our normalization 
is such that 

X, z) = 1. (41) 



Problemll(H): Show 



and use 



f 

Jo 



-a - Mixmi - Mix')) = -4 1 + in-^ 



1 o / _ _ X^ 

x^ 



to get (40) when orientations of Xqi and x^^ are averaged over. 
2.2.2 One soft gluon in the wavefunction 

The onium wavefunction having a single soft gluon, along with the heavy 
quark- antiquark, is calculated from the graphs in Fig. 4, where softness means 
Z2 << 1. In the large Nc limit the emission of a gluon changes a dipole into 
two dipoles. The original dipole, the heavy quark-antiquark pair, becomes 
two dipoles consisting of the heavy quark and the antiquark part of the gluon 
making up the first dipole and the quark part of the gluon along with the 
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heavy antiquark making up the second dipole. We illustrate the sum of the 
two graphs in Fig. 4 by the single graph of Fig. 5 where the dipole structure 
is emphasized. 





kj -k I-Z1-Z2 



(a) 



(b) 



FIG.4 

In order to derive the BFKL equation one needs only to calculate the 
graphs of Fig.4 in light-cone gauge where the polarization of the gluon is 



k 



(42) 



In the soft gluon approximation, which corresponds to the usual leading 
logarithmic approximation in which BFKL dynamics is formulated, one need 
only keep in (42) because the is a large quantity. Using the fact 

that the e_ polarization couples to the classical current of a high momentum 
quark or antiqark we can immediately write down the contribution of the 
graphs in Fig.4 to the wavefunction of the onium as 



^2/2/^2+ h+ 



{^^;^'\k, + k2, z) - ^^''\k^, z^)\ (43) 



where is the color matrix for a gluon of color a, and [l^/2k2+]'^ is the en- 
ergy denominator which is dominated by the soft gluon. Going to transverse 
coordinate space by 

^^'\xo2,X2l,z^,Z2) = J ^^|^e^^i-^oi+^^2-o.^(i)(A;,,fc2,^i,^2) (44) 
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X 



FIG.5 



one finds, using (43), 

TT V^02 2:21 J 

Squaring and taking a trace over colors gives 



2 ^2 



12-^02 



In obtaining (46) we have taken Cp = -/Vc/2 in the large Nc limit. 
Problem 12(E): Show 



E 



\Xq2 X2iy 



X, 



01 



™2 ™2 
•'^12-'^02 



Including a phase space factor 



one arrives at the factor 



$(0) 



(fx2 dZ2 



27r a;i2 2^02 
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as the probability of emitting a soft gluon, and where we have defined dy2 = 
dz2/z2- Eq.(47) is illustrated in Fig. 5. 



2.2.3 Arbitrciry numbers of soft gluons 

Once a single soft gluon emission has been calculated it is straightforward to 
add more, even softer gluons. The time sequence for soft gluons to appear in 
the wavefunction follows their longitudinal momentum. The hardest of the 
soft gluons appears first, then the next hardest of the soft gluons, etc. In each 
case the emission of a soft gluon creates a transition from a color dipole to two 
color dipoles with the probability factor given by (47). In order to express 
this most succinctly it is convenient to introduce a generating functional, Z, 
obeying 



Z{xoiY, u) = exp 



, , aNc 
u{xoi) + 



27r2 

J Ft XnoXio TT p 



TT p 

42X12 Jo n p 

R is a cutoff and indicates that one take X02 > P, X12 > p as the region of 
ingtegration. Eq.(48) is illustrated in Fig. 6. Z is the generating functional 
for soft gluons and has the following properties: 



FIG.6 



(49) 

Su{Axi) 5u{Ax2) 6u{Axn) 

gives, when multiplied by ^^'^^Xqi, Zi) the exclusive probability of having 
n — 1 soft gluons in the onium wavefunction with the soft gluons and the 
heavy quark- antiquark pair making up n dipoles of sizes Axi, Ax2, • • •. If one 
takes u = 1, rather than u = 0, the expression (49) becomes the inclusive 
probability of having n — 1 measured soft gluons, along with an arbitrary 
number of unmeasured gluons, in the onium wavefunction. Finally 
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Z {xoi, Y,u) =1 (50) 

u=l 

is probability conservation and is obtained by including the probability con- 
serving exponential factors in (48) which correspond to virtual loop correc- 
tions in the wavef unction ITJ 




FIG. 7 

Problem 13(H): Show that the virtual corrections, the terms in the expo- 
nentials in (48) can be obtained using probability conservation by integrating 
single gluon emission over the available (cutoff) phase space. That is 



Jo 71 P 



271^ Jr XQ2xi2 Jo 71 p 

In summary (48) represents the soft gluon wavefunction of a high energy 
heavy onium. The longitudinal momentum integrations have been done in 
a leading logarithmic approximation, while the transverse integrations have 
been handled exactly. The cutoff p will disappear when a physical problem 
is considered. The large Nc limit is esential to obtain (48). 

2.2.4 Scattering in the BFKL approximation 

In terms of the heavy onia light-cone wavefunctions, high energy onium- 
onium scattering is very simple and is reminiscent of partonic expressions in 
hard scattering. Suppose the scattering takes place in the center of mass. 
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Then the scattering proceeds by the interaction (scattering) of a single 
dipole in the left-moving onium with a single dipole in the right-moving 
onium as illustrated in Fig. 7. The equation is 



(51) 



where N{x, Y) is the number density of dipoles of size x, in a rapidity interval 
y, in the onium wavefunction. More explicitly 



N{x,Y) = j (fxoidzi^^^\xQi,zi)n{xQi,x,Y) 



with 



n{x,,,x,Y)^x^jdm 



5Z{xqi,Y,u) 



6u{x) 



u=l 



(52) 



(53) 



where d4>{x) indicates an integration over the orientations of the dipole di- 
rection X. n obeys the dipole version of the BFKL equation 



n{xoi, X, Y) — xS{x — Xoi)exp 



2aNc xoi 
Yin — 

P 



TT 



+ 



aNr. 



^2 ^2 



X [ dy 
Jo 

illustrated in Fig. 8 



exp 



.^-^iY-y)in^ 

TT P 



TT JR Xi2Xq2 

n{xi2:X,y) (54) 




FIG.8 

Problem 14 (M): Use (48) and (53) to derive (54). 

In order to solve (54) it is convenient to write it as an evolution equation 
in Y. Taking a Y-derivative on both sides of (54) gives 



d 2aN 

—n{xQi,x,Y)^ dx[2K{xQi,x[2)n{x[2,x,Y) (55) 

dY TT Jo 
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with 



1 /" 1 ijC 1 I 

K{xoi, = 75- / d'^X2S{xi2 - x'12) \ 2 °\ ~ 27r(5(x2 - XQ)£n{—) \ 

ZTT JR X'12 -^02 P ) 

(56) 

The hmit p — > can be taken in (56) leaving K a scale covariant kernel. 
Thus 

I rfa;i2i^(xoi, xi2)a;}f = x{i^)xlf''' (57) 
from the scale covariance. Exphcit calculation [6] gives 

Xiiy) = V'(l) - + - - ^^), (58) 

thus proving that K is the usual BFKL kernel in a different guise. Using 
(57) in (55) gives 

n(x„„x,y') = cr^(5i)'+-e^-W- (69) 

J —00 -^TT X 



Problem 15(E): Use n{xQi,x, 0) = xd{xQi — x) to how that c = 2 in (59). 



Probleml6(M): Use the fact that = is a saddle point of x(z^), along 
with the results x(0) = 2£n and x" = - 14C(3) to get 



n[Xoi,x,Y ) = — j====exp 



2x y|aiV,C(3)F 
Finally, using (52) and (60) in (51) gives 



nin 



2/xo 



UaNcC{3)Y 



a{Y) = IQna^R^ 



|aiV,C(3)F 



with 



and 



ap — 1 — 



AaN, 



en2 



TT 



(60) 



(61) 



(62) 
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R = lJ rfV^'c?2;iXoi$^°Ha;oi,2;i). (63) 

3 Phenomenology 

In this chapter we are going to try and confront the ideas of quark and 
gluon saturation with data. Before doing this we need to further refine 
the idea of saturation and relate it more precisely with that of unitarity. 
We begin our discussion by turning to soft hadronic interactions where the 
concept of parton does not apply but where there is significant evidence that 
cross sections have reached their unitary limit for small impact parameter 
scattering. 

3.1 Soft hadron-hadron scattering[14,15] 

Total and elastic cross sections are well represented in the Regge picture as 



(^tot 



ci{s/soY + C2/^fi/7o (64) 



with e ~ 0.1 and the constants Ci and C2 depending on the particular scat- 
tering being considered, e, however, appears to be universal. There is a 
major difficulty with a simple Regge pole picture, though, in that diffractive 
cross sections in proton-antiproton collisions grow much more slowly than 
expected ||16|| . The understanding of why diffractive cross sections grow more 
slowly at Tevatron energies than at lower energies is most easily seen in an 



impact parameter picture of high energy scattering 1 17]. One writes the to- 
tal, inelastic and elastic cross sections in terms of the S-matrix at an impact 
parameter b as 



atot = 2 J d%[l - S{h)\ 

= Jd%[l-S\b)] (65) 
aei = jd%[l-S{h)f. 
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FIG.9 

Equation (65) assumes that S{b) is real and it also neglects internal degrees of 
freedom in the proton relying on some average configuration in the proton as 
dominant in the scatterings described in (65). Using (65) to fit proton-proton 
and proton- antiproton scattering leads to a picture of S{b) as (schematically) 
illustrated in Fig.9. What is interesting in Fig. 9 is that the S-matrix is 
completely black, that is unitarity limits have been reached, for significant 
regions of b in high energy proton-antiproton collisions. In the regions which 
are completely black {S = 0) elastic and elastic cross sections are equal. 
Diffractive production can occur in the region where S goes from to 1 
but should not occur either when S is very near zero or when S is close to 
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one 



The increase in atot with energy occurs because some regions in b are 
changing from grey to black and because regions of larger b are going from 
white to grey. The region which contributes to diffraction grows less rapidly 
than those regions in b giving elastic and highly inelastic scattering. The 
Regge pole picture can be expected to be valid only when S is not too small. 
The regions where S is small correspond to multiple pomeron exchange in a 



Regge picture [19, 20 



3.2 Saturation of parton densities and unitarity 

3.2.1 Two useful reference frames in deep inealstic scattering 

Different frames of reference can be useful in understanding specific proper- 
ties of deep inelastic lepton-proton and lepton-nucleus scattering. For exam- 
ple, the Bjorken frame makes the view of the virtual photon as a probe of 
structure in the proton very natural. This picture of deep inelastic scattering 
is illustrated in Fig. 10 where p is very large at small x 




FIG. 10 



P~ (p + ^,o,o,p) 

q= (0,0,0,-Q). (66) 

The Bjorken variable x is equal to ^ in this frame which is close to the 
original Bjorken frame except that = so that the fact that the 7* is a 
local probe in x± at large Q is not manifest. However, in this frame the 
virtual photon can still be viewed as being absorbed by a quark over a short 
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period of time thus making manifest that F2 can be described in terms of the 
number of quarks in the proton. The frame just described, a Bjorken-hke 
or Breit-hke frame, is connected by a simple boost to another frame, which 
we shall call the dipole frame, where the process appears in time to proceed 
by the virtual photon breaking up into a quark- ant iquark pair which then 
scatters on the proton (or nucleus). In this frame, illustrated in Fig. 11, the 
momenta are 




FIG. 11 



TTl 

P~ + — ,0,0,p) 
2p 



g2-g2o,o,-g) 



(67) 



with X 



and where we choose q/Q greater than one but not 



too large. The idea is to have the virtual photon break up into a quark- 
antiquark pair but to keep all the structure and evolution in the proton 
rather than in a further evolution of the quark- ant iquark pair. This will be 
the case if q/Q is of moderate size. In this dipole frame the parton picture of 
deep inelastic scattering is not manifest, however, unitarity constraints are 
manifest pT|, p^ , p3| . In particular when a virtual photon breaks up into a 
quark- ant iquark pair having a definite transverse coordinate separation, x, 
and at a definite impact parameter, b, that pair scatters off the proton or 
nucleus with an S-matrix element S{b,x) which obeys \S{b,x)\ < 1. We shall 
see that S{b,x) = corresponds to saturation. Thus saturation (Bjorken 
frame) corresponds to blackness of the cross section (dipole frame). 
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3.2.2 Scattering off a nucleus (simple picture) 

To illustrate how quark saturation comes about consider deep inelastic scat- 
tering off a large nucleus in the dipole frame. We simplify the problem by only 
allowing the individual nucleons of the nucleus to interact with the quark- 
antiquark pair either through one gluon exchange, an inelastic reaction, or 
through two gluon exchange. We do not allow evolution in the individ- 
ual nucleons so one may view this model as an extension of the McLerran- 
Venugopalan model where now we have a (quantum) quark-antiquark pair 
interacting with the quasiclassical gluon field [^, |4|, |25|, |2^, ^ of the nu- 
cleus which we calculated in Chapter 1. However, since the nucleons in the 
nucleus do not have any evolution we can also simply take the nucleus at rest 
in which case the cross section takes the form of a virtual photon breaking 
into a quark-antiquark pair which then multiply scatters, via one and two 
gluon exchange, with the nucleons in the nucleus. This view is illustrated in 
Fig.l2. 



1 




and 



F2{x,Q')=Y.e}[xqf{x,Q')+xqf{x,Q')] (69) 
/ 

and the fact that aL/crr << 1 {ct and (Tl are the cross sections for trans- 
verse and longitudinal photons, respectively) one can directly write a a for- 
mula for the quark and antiquark distributions in terms of the cross section 
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for virtual photon scattering on the nucleus. In fact, we wish to go somewhat 
beyond this and write a formula for the measurement of a leading (current) 
jet in the deep inelastic scattering The equation isPSI 



with the photon wavefunction into a quark-antiquark pair of flavor / 



4a (21,^) = [^^^(l-z)[^^ + (l-^)^]gT/'e;^i(^?rtQV^(l-^))^, 

(71) 

and where the quark saturation momentum, Qg, obeys 

Ql = ^Ql (72) 

the various factors in (70) have a ready explanation. The factor of — 
comes from (68) relating a structure function (or parton distribution) to a 
virtual photon cross section. The wavefunctions give the amplitude, and com- 
plex conjugate amplitude, for a virtual photon to go into a quark-antiquark 
pair with momentum fractions z and 1 — z and at transverse coordinate 
separation Xi in the amplitude and X2 in the complex conjugate amplitude. 
The integration over (Pxi and (Px2 along with the e~*- '■-i~-2) factor fix the 
transverse momentum of the quark (or antiquark) to be i. The factors in the 
square bracket represent (respectively) no scattering of the pair, scattering 
in the amplitude and complex conjugate amplitude, scattering in the ampli- 
tude, and scattering in the complex conjugate amplitude. The exponential 
factors in the [ ] are recognizable as multiple scatterings if one identifies 

^ = 2VW^p^ (73) 

with p the usual nucleon density and cr the cross section of a quark-antiquark 
dipole of separation x on a nucleon. The | in front of the S in (70) is for nor- 
malization and means we are describing the jet rate for a single transversely 
polarized virtual photon. 
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If ^ is much less than Ql the integrals in (70) can be done giving 

dxqf Nc 



(74) 



which says there is ^ quarks, or antiquarks, per unit of transverse phase 
space in the wavefunction of the nucleus. This is quark saturation. This 
result follows from the unitarity of the dip ole- nucleus scattering. 
If one integrates (70) over i and b, the result 

xqf{x,Q'') + xqfix,Q^) = J d% J d\ dz z{l - z)\z^ + (1 - zf\ 



Kl[sjQ^x^z{l - z)]{l - e"^^'/^) (75) 
emerges. We have in fact reached what is essentially the Golec-Biernat and 



Wiisthoff model[^ a model of saturation successfully applied to both struc- 
ture functions and diffractive producton. The identification is complete if we 
identify Qs = l/Ro{x) and / d% = a^. Golec-Biernat and Wiisthoff do not do 
their analysis in impact parameter space but introduce a parameter setting 
the scale for cross sections, cxo, which is about 25 mb in their fits. They also 
parameterize Rl{x) = g^(^)'^ with A ~ 0.3 and xq ~ 3a;10~^. The fact 
that they are able to get such good fits of the data at low and moderate 

for F2 and for diffractive structure functions is, perhaps, the strongest 
evidence that saturation has been reached in the HERA energy range. 

Although (70) and (75) have been derived for scattering on a large nucleus 
the picture would seem to be pretty general. All the QCD dynamics is in the 
four factors in the [ ] in (70). The identification of the exponential with the 
S-matrix according to 

S{b,x) = e-^'^'/^ (76) 

means we may everywhere replace the gaussian exponential factors by the 
S-matrix and have a very general description of deep inelastic scattering, on 
either a proton or a nucleus, where unitarity is controlled. The essential 
physics which the Golec-Biernat and Wiisthoff model incorporates is 

S{h,x) = 1 when x^Ql/A « 1 (77) 
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S{b, x)=0 when x^Ql/A » 1 (78) 

for some scale Qg along with the statement that the b dependence is on the 
scale of a fermi for protons and several fermis for nuclei. 

3.2.3 Is there evidence for saturation? 

Although the Golec-Biernat and Wiisthoff model gives the strongest evidence 
for saturation I would like to go over the argument in a very qualitative way 
because the issue is so important and because the physics argument, though 
not technically difficult, is subtle. In general fits to F2{x, Q^) using a second 
order DGLAP formalism work very well in the HERA regime. However, 
there are signs that all is not well with such fits in the region of very small x 
and moderate values of Q^. In addition there is evidence that something like 
saturation is needed in order to be consistent with the energy dependence of 
diffractive structure functions measured at HERA. 

We begin with F2{x, Q"^). Sometime ago A. Caldwell pO|, [31[] presented ev- 



idence for a change in the Q^— dependence of F2 at very small x and mod- 
erate . The HERA data is shown in Fig. 13 where an unexpected turnover 
in g^^2 occurs at moderate and small x. Data through a similar region of 
Q^, but at much larger values of x from fixed target experiments are shown 
in Fig. 14 where there is no turnover visible at moderate Q^. This turnover 
is well described by a formula like that given in (75), that is by the model 
of Golec-Biernat and Wiisthoff. Let us here try to describe what may be 



happening in very simple terms [32, 33 . 
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FIG. 13 

In the dipole frame the virtual photon breaks up into a quark-antiquark 
pair of relative momentum 2k which then scatters, inelastically, on the proton 
as illustrated in Fig. 15. We may write 



where 



dF2 oc Q^dPr{P)a,4k) 



dm') oc ^ 



(79) 



(80) 



is the probability that the virtual photon break up into the quark-antiquark 
pair in the momentum range c/fc^. crj„ is the cross section for the pair to 
interact with the photon. In perturbation theory 



a 



C7f-*(fc)(XpXG(x,i^) 



(81) 
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FIG. 14 



so that cancels out in (79) and one is left with the usual d^k^ integral 
giving scahng violations. In this case Q^f^ axG{x,Q'^) and one would 
interpret the turnover in the Caldwell plot as a (surprising) decrease in the 
gluon distribution at very low values of x. However, if is in the saturation 
regime 

at:\k) = TiRl (82) 

with Rq about the proton's radius, since the unitarity limit is being reached. 
This means that Q^f^ ~ Q^Rq and the turnover at small is seen as a 
natural result of unitarity, or saturation. It would be expected that unitarity, 
and saturation, should occur earlier for small impact parameter interactions 
than for large impact parameter interactions. It would be very interesting to 
be able to do analysis of the Q^— dependence of structure functions at very 
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FIG. 15 



small function of the impact parameter of the collision. 

The other process for which there is some evidence for saturation is in the 
diffractive cross section at HERA [p9|, ^ ^ The picture again is as in 



Fig. 15 but where now the proton, p, must exchange at least two gluons, and 
remain in a color singlet, in order that a large rapidity gap emerge with the 
proton typically appearing in the final state at a small momentum transfer. 
One can write 

daoif oc dPr{k)aeiik) (83) 
with dPr as in (80). In perturbation theory 

aTcx[^,xG{x,e)r (84) 

with X ~ M'^/s and with M the mass of the diffractive state coming from 
the quark-antiquark pair. Using (84) in (83) leads to a quadratically di- 
verging infrared integration d]^ /{J^Y, and this suggests that the process of 
diffraction at HERA is a soft process. This was the common wisdom among 
theorists until a few years ago. What could happen to make the process 
semihard? If the saturation momentum is as large as a GeV or a few GeV 
then (84) applies for transverse momenta above the saturation regime. If \k\ 
is below the saturation regime one should use 

<f oc T^Rl (85) 
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in (83). This means that the saturation momentum serves as the infrared 
cutoff for the "divergent" integral when (82) is used in (81). If the saturation 
momentum is in the semihard regime (1-2) GeV then the resulting process 
is not a soft process and one might expect an energy dependence which is 
stronger than that given by the soft pomeron. In fact the data support a 
fairly strong energy growth. If one writes 

oc x--* (86) 

ax 

then both ZEUS and HI find 7 to be about a factor of 2 larger than that 
expected from soft pomeron exchange. What is unique about the scattering 
of a virtual photon on a hadron is that the virtual photon wavefunction is 
dominated by high transverse momentum states, in contrast to hadrons. If 
the unitarity limit, saturation, is reached in some range of transverse momen- 
tum the process simply proceeds to a higher region of transverse momentum. 
In diffraction of virtual photons the most important region of transverse 
momenta in the virtual photon's wavefunction should be just above the sat- 
uration momentum. 
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